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ABSTRACT 

We introduce a novel method for weak-lensing measurements, which is based on a mathemat- 
ically exact deconvolution of the moments of the apparent brightness distribution of galaxies 
from the telescope's PSF. No assumptions on the shape of the galaxy or the PSF are made. The 
(de)convolution equations are exact for unweighted moments only, while in practice a com- 
pact weight function needs to be applied to the noisy images to ensure that the moment mea- 
surement yields significant results. We employ a Gaussian weight function, whose centroid 
and ellipticity are iteratively adjusted to match the corresponding quantities of the source. The 
change of the moments caused by the application of the weight function can then be corrected 
by considering higher-order weighted moments of the same source. Because of the form of 
the deconvolution equations, even an incomplete weighting correction leads to an excellent 
shear estimation if galaxies and PSF are measured with a weight function of identical size. 

We demonstrate the accuracy and capabilities of this new method in the context of weak 
gravitational lensing measurements with a set of specialized tests and show its competitive 
performance on the GREAT08 challenge data. A complete C++ implementation of the method 
can be requested from the authors. 
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1 INTRODUCTION 

Shear estimation from noisy galaxy images is a challenging task, 
even more so for the stringent accuracy requirements of upcoming 
cosmic shear surveys. Existing methods can achieve multiplicative 
errors in the percent range (Bridle et al. 2010), but to exploit the sta- 
tistical power of the next-generation surveys errors in the permille 
range or even below are requested (Amara & Refregier 2008). 

Shear estimates can be achieved in a model-based or in a 
model-independent fashion. For instance, Lensfit (Miller et al. 
2007) compares sheared and convolved bulge-disk profiles to the 
given galaxies. Model-based approaches often perform excellently 
for strongly degraded data because certain implicit or explicit priors 
keep the results within reasonable bounds, e.g. the source ellipticity 
smaller than unity. On the other hand, when imposing these priors 
to data, whose characteristics differ from the expectation, these ap- 
proaches may also bias the outcome. 

Model-independent approaches do not - or at least not as 
strongly - assume particular knowledge of the data to be analyzed. 
They should therefore generalize better in applications, where pri- 
ors are not obvious, e.g. on the intrinsic shape of lensed galaxies. 
The traditional KSB method (Kaiser et al. 1995) forms a shear es- 
timator from the second-order moments of lensed galaxy images. 
When doing so, it is not guaranteed that reasonable shear estimates 
can be achieved for each galaxy. Consequently, KSB requires a 
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careful setup, which is adjusted to the characteristics of the data 
to be analyzed. KSB furthermore employs strong assumptions on 
the PSF shape, which are not necessarily fulfilled for a given tele- 
scope or observation (Kuijken 1999). As we have shown recently, 
KSB relies on several other assumptions concerning the relation 
between convolved and unconvolved ellipticity as well as the rela- 
tion between ellipticity and shear, neither of which hold in practice 
(Viola et al. 2011). 

In this work we present a novel method for shear estimation, 
which maintains the strengths of model-independent approaches 
by working with multipole moments, but does not suffer from the 
KSB-shortcomings mentioned above. 



2 THE Deimos METHOD 

The effect of gravitational lensing on the surface brightness dis- 
tribution G(x) of a distant background galaxies is most naturally 
described in terms of the moments of the brightness distribution, 

{G}yS J rf 2 xG(x)x',4 (1) 

for which we introduce a tensor-like notation. The effect of the re- 
duced shear g is contained in the change of the complex ellipticity 

{Gh,o - {G)o, 2 + 2i{G} u 
X " {G) 2 , + (G)o, 2 
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with respect to its value before lensing, 

X ~ 2g + g 2 X* 



1 + |g| 2 - 2%{g X *) 



(3) 



(e.g. Bartelmann & Schneider 2001). Unfortunately, we do not 
know x\ which would allow solving directly for g given x- 
Furthermore, (at least) two observational complications alter the 
source's moments, typically much more drastically than lensing: 
convolution with the PSF and any means of noise reduction - nor- 
mally weighting with a compact function - to yield significant mo- 
ment measurements. Hence, \ is not directly accessible and needs 
to be estimated by properly accounting for the observation effects. 
We leave the treatment of weighting for section 3 and start with the 
derivation of the change of the moments under convolution. 

Any square-integrable one-dimensional function G(x), has an 
exact representation in Fourier space, 



G(x) -> tf> G (k) : 



dxG{x)e k 



(4) 



In the field of statistics, ipc is often called the characteristic func- 
tion of G and has a notable alternative form 1 



1> G (k) = Y\G} n — 



(iky 
i 



(5) 



which provides a link between the Fourier-transform of G and its 
moments (G)„, the one-dimensional pendants to Equation 1. We 
can now employ the convolution theorem, which allows us to re- 
place the convolution by a product in Fourier-space, i.e. by a prod- 
uct of characteristic functions of G and of the PSF kernel P, 



G*(x) 



= f dx G(x 



)P(x-x') 



(6) 



For convenience, we assume throughout this work the PSF to be 
flux-normalized, (P)o = 1. Considering Equation 5, we get 

n n 

(ik)"' , (ik)' 



y. y.{G) m — — [p} n - m - — — 

n m 



(7) 



i QkT 

I n\ ' 



where we applied the Cauchy product in the second step. The ex- 
pression in square brackets on the last line is by definition the de- 
sired moment, 



{G*}„ =^(:)fGUP}„- 



(8) 



Hence, we can now express a convolution of the function G with 
the kernel P entirely in moment space. Moreover, even though the 
series in Equation 5 is infinite, the order of the moments occurring 
in the computation of (G*)„ is bound by n. This means, for calculat- 
ing all moments of G* up to order n, the knowledge of the same set 
of moments of P and G is completely sufficient. This results hold 
for any shape of G and P as long as their moments do not diverge. 
For non-pathological distributions, this requirement does not pose 
a significant limitation. 

An identical derivation can be performed for two-dimensional 

1 The summation indices in this work all start with zero unless explicitly 
noted otherwise. 



Table 1. Equations for deconvolving all moments up to order n = 2. The 
shown equations are specializations of Equation 9. 
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moments, which yields the change of moments of G(x) under con- 
volution with the kernel P(x) (Flusser & Suk 1998): 



(9) 



Deconvolution 

To obtain the deconvolved moments required for the shear estima- 
tion via the ellipticity x> we need to measure the moments up to 
second order of the convolved galaxy shape and of the PSF kernel 
shape. Then we can make use of a remarkable feature of Equation 9, 
which is already apparent from its form: The impact of convolu- 
tion on a moment of order i + j = n is only a function of uncon- 
volved moments of lower order and PSF moments of at most the 
same order. We can therefore start in zeroth order, the flux, which 
only needs to be corrected if the PSF is not flux-normalized. With 
the accurate value of the zeroth order and the first moments of the 
PSF, we can correct the first-order moments of the galaxy, and so 
on (cf.Table 1). The hierarchical build-up of the deconvolved mo- 
ments is the heart of the Deimos method (short for deconvolution in 
moment space). 

It is important to note and will turn out to be crucial for weak- 
lensing applications that with this deconvolution scheme we do not 
need to explicitly address the pixel noise, which hampers most 
other deconvolution approaches in the frequency domain, simply 
because we restrict ourselves to inferring the most robust low-order 
moments only. 



3 NOISE AND WEIGHTING 

In practice, the moments are measured from noisy image data, 

I(x) = G(x) + N(x), (10) 

where the noise N can be considered to be independently drawn 
from a Gaussian distribution with variance a\, i.e. (N(Xt)N(Xj)) = 
o~ldjj for any two positions x, and Xj. According to Equation 1, the 
image values at large distances from the galactic center have the 
largest impact on the (/)„ if n > 0. For finite and compact brightness 
distributions G, these values are dominated by the noise process in- 
stead of the galaxy, whose moments we seek to measure. Conse- 
quently, centered weight functions W of finite width are typically 
introduced to limit the integration range in Equation 1 to regions in 
which / is mostly determined by G, 



Ux) = W(x)I(x). 



(11) 
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Table 2. Correction terms for deweighting moments of order n = i + j. The 
deweighted moments {Idw}i,i are given by the sum of the correction terms 
up to the limiting order n w . 

n w correction terms 
{I w )ij 

2 JJJ [ci{/„,),- +2 J - 4 e 2 |/,r),+ I J+l + C2{/ w }ij +2 ] 

4 ^4 fcf{/ u -}/+4j - 8cie{/ w } i+3 j + i + [2cic 2 + 16ef],{/w}>+2j+2 - 

6 ^6[ei{Wi+6,j - 12c 2 e 2 {/ u .), +5J+ i + [3f2 C2 + 48cie|] j/„,) i+4J+2 - 
[24c lC2 e 2 + 64e 2 3 ]{/„,|, +3 j + 3 + [3 c,c\ + 48 c 2 ej] {I w \ 
12c|e 2 {/ w |, + ij +5 + (^{/wJij+s] 



A typical choice for W is a circular Gaussian centered at the galac- 
tic centroid, 

W(x) = exp(-^). (12) 

Alternatively, one can choose to optimize the weight function to 
the shape of the source to be measured. Bernstein & Jarvis (2002, 
see their section 3.1.2) proposed the usage of a Gaussian, whose 
centroid x c , size s, and ellipticity e are matched to the source, such 
that the argument of the exponential in Equation 12 is modified 
according to 

x - x '=( 1 -;r i?J (x -^ (i3) 

As such a weight function represents a matched spatial filter, it op- 
timizes the significance and accuracy of the measurement if its pa- 
rameters are close to their true values. This can, however, not be 
guaranteed in presence of pixel noise, but we found the iterative al- 
gorithm proposed by Bernstein & Jarvis (2002) to converge well in 
practice and therefore employ it to set the weight function within 
the Deimos method. 

Unfortunately, a product in real space like the one in Equa- 
tion 1 1 translates into a convolution in Fourier-space. We therefore 
have to expect some amount of mixing of the moments of /„,. Even 
worse, an attempt to relate the moments of /„ to those of / leads to 
diverging integrals. Hence, there is no exact way of incorporating 
spatial weighting to the moment approach outlined above. On the 
other hand, we can invert Equation 11 for / = I w /W and expand 
1/ IV in a Taylor series around the center at x = 0, 

2 

W~\x) * W-'(0) - 1V'(0)[£ c k x\ + 4e 2 x,x 2 ] + 

k=l 

1 2 

- W"(0)[2 CkC,x 2 k xf - 8e 2 ^ c k x 2 kXl x 2 + (4e 2 x,x 2 ) 2 ], 

k,l=i k=l 

(14) 

where we employed W'(x) = and ci, 2 = (1 + ei) 2 + e|. We 
introduce the parameter n„, as the maximum order of the Taylor 
expansion, here n w = 4. Inserting this expansion in Equation 1, 
we are able to approximate the moments of / by their deweighted 
counterparts {I dw } . For convenience we give the correction terms for 
orders n u , < 6 in Table 2. This linear expansion allows us to correct 
for the weighting-induced change in the moments of a certain order 
n by considering the impact of the weight function on weighted 
moments up to order n + n w . 



3.1 Deweighting bias 

The truncation of the Taylor expansion constitutes the first and only 
source of bias in the Deimos method. The direction of the bias is ev- 
ident: As the weight function suppresses contributions to the mo- 
ments from pixel far away from the centroid, its employment re- 
duces the power in any moment by an amount, which depends on 
the shape - particularly on the radial profile - of the source and the 
width s. Additionally, if the ellipticity e was misestimated during 
the matching of W, the measured ellipticity of the source x be- 
fore and after deweighting will be biased towards e. In the realistic 
case of noisy images, which we address in more detail in subsec- 
tion 3.2, e can be wrong in two ways: statistically and systemat- 
ically. The centered Gaussian distribution of the pixel noise leads 
to a centered Cauchy-type error distribution of both components of 
e, i.e. the statistical errors of e and therefore also;^ have vanishing 
mean. The systematic errors stem from the application of a compact 
weight function to measure e, which for any finite s constitutes a 
removal of information. This leads to deviations of the measured 
e from the true e if e.g. the centroid is not determined accurately 
or the ellipticity changes with radius. For individual objects, these 
deviations are impossible to quantify precisely - as this would re- 
quire knowledge of the true observed morphology - and hamper 
all weak-lensing measurements unless an appropriate treatment is 
devised (e.g. Hosseini & Bethge 2009; Bernstein 2010). For large 
ensembles of galaxies, effective levels of systematic errors can be 
assessed by analyzing dedicated simulations (cf. section 4). 

We investigate now the DEiMos-specific systematic impact of a 
finite n w on the recovery of the deweighted moments. For the exper- 
iments in this section we simulated simple galaxy models following 
the Sersic profile 

p.Wocexpj-fo,,,^)"" 1 -!]}, (15) 

where n s denotes the Sersic index, and^ f the effective radius, while 
the PSFs are modeled from the Moffat profile 

Pm (r)K(l+ar 2 )-P, (16) 

where a = (2 I//? - l)/(FWHM/2) 2 sets the width of the profile and 
p its slope. Both model types acquire their ellipticity according to 
Equation 13 

In the top panel of Figure 1 we show the error after deweight- 
ing a convolved galaxy image from a matched elliptical weight 
function as a function of its size s. As noted above, the bias is al- 
ways negative and is clearly more prominent for the larger disk- 
type galaxy (circle markers). As the Taylor expansion becomes 
more accurate for n„. — > <x> or s — > oo, the bias of any moment 
decreases accordingly. 

An important consequence of the employment of a weight 
function with matched ellipticity is that the bias after deweighting 
does only very weakly depend on the apparent ellipticity, i.e. all 
moments of the same order are biased by the same relative factor 
A(n, s). This means any ratio of such moments remains unbiased. 
This does not guarantee that the ellipticity is still unbiased after 
the moments have passed the deconvolution step, which is exact 
only for unweighted moments. On the other hand, the particular 
form of the equations in Table 1 becomes important here: If we as- 
sume well-centered images of the galaxy and the PSF and a negli- 
gible error of the source flux (G)o,o> none of which is guaranteed for 
faint objects, the deconvolution equations for the relevant second- 
order moments only mix second-order moments. If furthermore 
A G (2, s) = A/=(2, s), the ellipticity x ( c f- Equation 2) will remain 
unbiased after deconvolution even though the moments themselves 
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The variance of the weighted moments is given by 



Figure 1. Weighting-induced bias. Relative error of the PSF-convolved and 
deweighted moment [Idw}l,l (top) and of the estimated ellipticity after de- 
convolution (bottom) as a function of weight function size s. Colors indi- 
cate the correction order n w , while markers denote the galaxy model (cir- 
cles for model 1, diamonds for model 2). The galaxy models are simulated 
as Sersic profiles with the following parameters: e s = (0.1,0.4), n, = 
1(4), R e = 3(1.5) for model 1 (2). The PSF was of Moffat-type with 
e = (0.05,0.05), /3 = 3, FWHM = 3. The weight functions of galaxy 
and PSF had the same size s, but individually matched ellipticities. 



were biased. The aforementioned condition holds if the radial pro- 
files of PSF and galaxy are similar within the weight function, in 
other words, if the galaxy is small. This behavior can clearly be 
seen in the bottom panel of Figure 1 , where the ellipticity estimate 
of the smaller elliptical galaxy (diamond markers) has sub-percent 
bias for n w > 2 and s > 3. The estimates for the larger galaxy 
are slightly higher because |Ag(2, s)| > |A/>(2, s)\, i.e. the deconvo- 
lution procedure overcompensates the PSF-induced change of the 
moments. However, sub-percent bias is achieved for n w > 4 and 
s > 5. 

For large galaxies, it might be advantageous to adjust the sizes 
of galaxy and PSF independently as this would render A G (2, s G ) 
more comparable to A P (2, s P ). However, we found employing a 
common size s for both objects to be more stable for the small 
and noisy galaxy images typically encountered in weak-lensing ap- 
plications. We therefore adjust the size s such as to allow an opti- 
mal measurement of the deweighted PSF moments. Since the main 
purpose of the weighting is the reduction of noise in the measured 
moments, one could improve the presented scheme by increasing 
s for galaxies with larger surface brightness such as to reduce the 
bias when the data quality permits. 



3.2 Deweighting variance 

Being unbiased in a noise-free situation does not suffice for a prac- 
tical weak-lensing application as the image quality is strongly de- 
graded by pixel noise. We therefore investigate now the noise prop- 
erties of the deweighted and deconvolved moments. 



r; J dx W 2 (x)xfx 



(17) 



since the noise is uncorrected and has a vanishing mean. It is ev- 
ident from Table 2 that the variance of the deweighted moments 
increases with the number of contributing terms, i.e. with n w . Less 
obvious is the response under changes of s. While each moment ac- 
cumulates more noise with a wider weight function, the prefactors 
of the deweighting correction terms is proportional to s~" w such that 
their impact is reduced for larger s. 

To quantitatively understand the impact of n w and s in a fairly 
realistic scenario we simulated 10,000 images of the galaxy models 
1 and 2 from the last section. We drew their intrinsic ellipticities 
from a Rayleigh distribution with <T\ c s\ = 0.3. Their flux was fixed 
at unity, and the images were degraded by Gaussian pixel noise 
with variance cr\. We ran Deimos on each of these image sets with 
a fixed scale The results are presented in Figure 2, where we 
show the dispersion of the measured x m units of the dispersion 
of X s - From the left panel it becomes evident that the attempt of 
measuring unbiased ellipticities (large n w or s) comes at the price 
of increased noise in the estimates. Considering also Figure 1, we 
infer that in this bias-variance trade-off small values of s and large 
values of n„. should be favored since this provides estimates with 
high accuracy and a moderate amount of noise. 

In the right panel of Figure 2 we show the estimator noise as 
function of the pixel noise. Equation 17 suggests that there should 
be a linear relation between these two quantities, which is roughly 
confirmed by the plot. Additional uncertainties in the moment mea- 
surement - caused by e.g. improper centroiding - and the non- 
linear combinations of second-order moments to form x lift ine 
actual estimator uncertainty beyond the linear prediction. 

Even though the true errors of x mav not exactly follow the 
linear theory, we will now exploit the fairly linear behavior to form 
error estimates. We can express the deweighting procedure as a ma- 
trix mapping, 



U,/ 



D • {!„.) 



(18) 



where {I w } and (I Av ) denote the vectors of all weighted and 
deweighted moments, and D encodes the correction terms of Ta- 
ble 2. The diagonal covariance matrix S K . of the weighted moment 
variances given by Equation 17 is related to the covariance matrix 
of the deweighted moments by 

S rf = D S„ D T , (19) 

from which we can obtain the marginalized errors by 

(t7 2 ({/ rfw } ij .))" 1 =(S^ 1 ) ifcjt) (20) 

where k denotes the position of the moment {/Jy in the vector 
. Under the assumptions mentioned above, also the deconvo- 
lution can be considered a linear operation, at least up to order 
2, so that we can extend the error propagation even beyond this 
step: If we neglect errors in the PSF moments, the errors of the 
deconvolved moments (up to order 2) are identical to those of the 
deweighted ones. We can therefore estimate the errors of all quan- 
tities based on deconvolved moments directly from Equation 20. 



4 SHEAR ACCURACY TESTS 

So far, we were concerned with the estimation of ellipticity. To test 
the ability of our new method to estimate the shear, we make use 
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Figure 2. Noise of the deweighted and deconvolved ellipticity estimate as a 
function of the weight function size s (left) and of the standard deviation of 
the pixel noise <r„ (right). The color and marker code is explained in Fig- 
ure 1. The pixel noise is given in units of 10~ 3 for flux-normalized sources. 
cr„ = 8 is close to the detection limit for this source model. 



of the reference simulations with realistic noise levels from the 
GREAT08 challenge (Bridle et al. 2010). As the shear values in 
these simulations are fairly low, we employ the linearized version 
of Equation 3, corrected by the shear responsivity of the source en- 
semble, 



(X) 



(21) 



(e.g. Massey et al. 2007), without any further weighting of indi- 
vidual galaxies, to translate Demos ellipticity measures into shear 
estimates. The dispersion cr? is measured from the lensed and noisy 
galaxy images and hence only coarsely describes the intrinsic shape 
dispersion (cf. Figure 2). We are aware of this limitation and ver- 
ified with additional simulations that it introduces sub-percent bi- 
ases for the range of shears and pixel noise levels we expect from 
the GREAT08 images. 

We inferred the weight function size s = 4 and the correction 
order ;i„. = 4 from the optimal outcome for a set with known shears. 
The actual GREAT08 challenge data comprises 9 different image 
sets, which differ in the shape of the PSF, the signal-to-noise ratio, 
the size, and the model-type of galaxies. For each of these branches, 
there are 300 images with different values of shear. We performed 
the Deimos analysis of all images keeping the weighting parameters 
fixed to the values inferred before. The results are shown in Figure 3 
in terms of the GREAT08 quality metric Q (see eqs. 1& 2 in Bridle 
et al. 2010) and of the multiplicative shear accuracy parameters ra, 
obtained from a linear fit of the shear estimates g, to the true shear 
values g, (Heymans et al. 2006; Massey et al. 2007), 



: m, gj + c, . 



(22) 



From Figure 3 we clearly see the highly competitive perfor- 
mance of Deimos with a typical Q > 200 in all but two branches. 
Single-component galaxy models yield a particularly large Q- 
value, probably because the bulge-only models are the most com- 
pact ones and thus favor the setting of a constant s for PSF and 
galaxies. In terms of Q, there is no change between the centered 
and the off-centered double-component galaxy models, but both ra, 
drop for the off-centered ones. As such galaxy shapes have vari- 
able ellipticity with radius and Deimos measures them with a fixed 
weight function size, we interpret this as a small but noticeable 
ellipticity-gradient bias (Bernstein 2010). 
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Figure 3. GREAT08 Q metric and multiplicative shear accuracy m, for the 
nine different branches of the GREAT08 challenge data with realistic noise 
levels. In each panel, the scale on the left describes the values of Q and 
the scale on the right the values of m, . The dotted line denotes m, = 0. 
The dashed lines show the Q- value before we adjusted the weight-function 
matching and deweighting parameters to the source characteristics of the 
branch. 



The response to changes in the PSF shape is a bit more wor- 
risome and requires explanation. The fiducial PSF had e\ > and 
the opposite is true for the rotated one. From all panels of Figure 3 
we can see that typically \mi\ < \m 2 \. Such a behavior has already 
been noted by Massey et al. (2007): Because a square pixel appears 
larger in diagonal direction than along the pixel edges, the moment 
{/)n and hence e 2 suffer more strongly from the finite size of pixels. 
From the discussion in subsection 3.1, we expect a certain amount 
of PSF-overcompensation for small weighting function sizes. As 
the PSF shape is most strongly affected by pixelation, the overcom- 
pensation boosts preferentially those galaxy moments, which align 
with the semi-minor axis of the PSF. In general, a larger PSF - or a 
larger PSF ellipticity - improves the shear estimates. It is important 
to note, that, as in all other panels, the residual additive term c,- was 
negligible for all PSF models. 

The response to changes in S/N or galaxy size is more dra- 
matic: Particularly the branches 7 (low S/N) and 9 (small galaxies) 
suffer from a considerable shear underestimation. This is not sur- 
prising as also most methods from Bridle et al. (2010) showed their 
poorest performance in these two sets. Since the Q metric strongly 
penalizes poor performance in single GREAT08 branches, the over- 
all Q = 7.7 for this initial analysis. 

As this is the first application of Deimos to a weak-lensing test 
case, we allowed ourselves to continue in a non-blind fashion in or- 
der to work out how the Deimos estimates could be improved. Ap- 
parently, problems arise when the galaxies are small or faint. The 
obvious solution is to shrink the weight function size. As discussed 
in section subsection 3.2, improper centroiding plays an increasing 
role in deteriorating shear estimates for fainter galaxies. We there- 
fore split the weight-function matching into two parts: centroid de- 
terminations with a small weight function of size s c , and ellipticity 
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determination with s > s c . By choosing s c = 1.5 and s = 2.5, 
we could strongly improve the performance for branches 7 and 9. 
Given the high S/N of branch 6, we decided to rerun these images 
with n w = 6, which yielded another considerable improvement. 
With these modifications to the weight-function matching and the 
deweighting parameters, Deimos estimates achieve Q = 112, simi- 
larly to Lensfit with Q = 119, at a fraction of the runtime (0.015 
seconds per GREAT08 galaxy). We emphasize that this is a some- 
what skewed comparison as we had full knowledge of the simula- 
tion characteristics. However, the changes to the initial analysis are 
modest and straightforward. In particular, they depend on galactic 
size and magnitude only, and not on the true shears. 

Given the bias-variance trade-off from the deweighting pro- 
cedure, the outcome of this section also clearly indicates that a 
simple one size fits all approach is not sufficient to obtain highly 
accurate shear estimates from Deimos. For a practical application, 
a scheme to decide on n w , s c , and s for each galaxy needs to be in- 
corporated. Such a scheme can easily be learned from a small set of 
dedicated simulations, foremost because the Deimos results depend 
only weakly on PSF and intrinsic galaxy shape. 



5 COMPARISON TO OTHER METHODS 

Because of the measurement of image moments subject to a 
weighting function, Deimos shares basic ideas and the computa- 
tional performance with the traditional KSB-approach (Kaiser et al. 
1995). In contrast to it, Deimos does not attempt to estimate the 
shear based on the ellipticity of single galaxies 2 , nor does it need 
to assume that the PSF can be decomposed into an isotropic and 
an anisotropic part, which introduces residual systematics into the 
shear estimation if the anisotropic part is not small (Kuijken 1999). 
Deimos rather offers a mathematically exact way of deconvolving 
the galaxy moments from any PSF, thereby circumventing the prob- 
lems known to affect KSB (see Viola et al. (201 1) for a recent dis- 
cussion of the KSB shortcomings). Its only source of bias stems 
from the inevitably approximate treatment of the weight function, 
which requires the measurement of higher-order image moments. 
Since Deimos measures all moments with the same weight func- 
tion (instead of with increasingly narrower higher derivatives of 
the weight function), these higher-order correction terms suffer less 
from pixelation than those applied in KSB. However, as we could 
see in section 4, pixelation affects the Deimos measurements, and 
an analytic treatment of it is not obvious. 

The treatment of the convolution with the PSF on the basis of 
moments is very close to the one known from shapelets (Refregier 

6 Bacon 2003; Melchior et al. 2009). However, Deimos does not 
require the time-consuming modeling process of galaxy and PSF, 
and hence is not subject to problems related with insufficient mod- 
eling of sources, whose apparent shape is not well matched by a 
shapelet model of finite complexity (Melchior et al. 2010). 

In the RRG method (Rhodes et al. 2000), the effect of the 
PSF convolution is also treated in moment space. Furthermore, an 
approximate relation between weighted and unweighted moments 
is employed, which renders this approach very similar to the one 
of Deimos. The former differs in the employment of the KSB-like 
anisotropy decomposition of the PSF shape. 

As mentioned in section 3, Deimos makes use of the same it- 
erative algorithm as ELLIPTO (Bernstein & Jarvis 2002) to define 

2 This demands setting x s = i n the non-linear Equation 3, which is only 
true on average but not individually. 



the centroid and ellipticity of the weight function. The latter addi- 
tionally removes any PSF anisotropy by applying another convolu- 
tion to render the stellar shapes circular, which is not necessary for 
Deimos. 

The recently proposed FDNT method (Bernstein 2010) decon- 
volves the galaxy shape from the PSF in the Fourier domain, and 
then adjusts centroid and ellipticity of the coordinate frame such 
that the first-order moments and the components of the ellipticity 
- formed from second-order moments - vanish in the new frame. 
FDNT restricts the frequencies considered during the moment mea- 
surement to the regime, which is not suppressed by PSF convolu- 
tion. Because of the shearing of the coordinate frame, additional 
frequencies need to be excluded, whereby the allowed frequency 
regime further shrinks. This leads to reduced significance of the 
shear estimates for galaxies with larger ellipticities. Furthermore, 
FDNT requires complete knowledge of the PSF shape. In con- 
trast, Deimos does not need to filter the data, it extracts the lensing- 
relevant information from the low-order moments of the galaxy and 
PSF instead. These differing aspects indicate that Deimos should be 
more robust against pixel noise. It should also be possible to in- 
corporate the correction for ellipticity-gradient bias suggested by 
Bernstein (2010) in the Deimos method. 



6 CONCLUSIONS 

For the presented work, we considered the most natural way of 
describing the effects of gravitational lensing to be given by the 
change of the multipole moments of background galaxies. We di- 
rectly estimate the lensed moments from the measured moments, 
which are affected by PSF convolution and the application of a 
weighting function. For the PSF convolution we derive an analytic 
relation between the convolved and the unconvolved moments, 
which allows an exact deconvolution and requires only the knowl- 
edge of PSF moments of the same order as the galaxy moments to 
be corrected. The weighting-induced changes of moments cannot 
be described analytically, but for smooth weight functions a Taylor 
expansion yields approximate correction terms involving higher- 
order moments. 

We showed that the residual bias of the deweighted moments 
stemming from an incomplete weighting correction is modest. 
Moreover, choosing a weight function with matched ellipticities but 
same size for measuring stellar and galactic moments yields ellip- 
ticity estimates with very small bias even for rather small weight- 
ing function sizes, which are required to reduce the impact of pixel 
noise to a tolerable level. In this bias-variance trade-off, Deimos 
normally performs best with high correction orders n w at small 
sizes j, but data with high significance may need a different setup. 
The choice of these two parameters is the trickiest task for a Deimos 
application, but can be easily addressed with a dedicated simula- 
tion, which should resemble the size and brightness distribution of 
sources to be expected in the actual data. Other properties of the 
sources, like their ellipticity distribution or, more generally, their 
intrinsic morphology, do not need to be considered as the measure- 
ment of moments does neither imply nor require the knowledge of 
the true source model. 

There are certain restrictions of the method to bear in mind: 

(i) Setting s to be the same for galaxies and the PSF works best 
for small galaxies, whose shape is dominated by the PSF shape. 

(ii) Changes of the shape at large radii would fall outside of the 
weight function and hence be ignored. When present in the PSF 
shape, this could lead to a residual PSF contamination, but can be 
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cured by increasing the scale of the weight function at the expense 
of larger noise in the galaxy moments. When present in galactic 
shapes, the results become susceptible to ellipticity-gradient bias. 

(iii) Direct measurement of the moments from the pixel values is 
inevitably affected by pixelation. For small, potentially undersam- 
pled shapes this leads to biased moment and ellipticity measures 
and acts more strongly in diagonal direction, i.e. on e 2 . 

(iv) The noise on the ellipticity estimates based on image mo- 
ments is not Gaussian, nor does it propagate easily into the shear es- 
timate. When dominant, it can create substantial biases of its own. 

Only the first of these restrictions exclusively applies to Demos, 
the others are present in all non-parametric methods, which work 
directly on the pixelated image. Model-based approaches could 
replace the coarsely sampled moment measurements by ones ob- 
tained from the smooth models. 

Further work is required to choose the deweighting param- 
eters, to account for pixelation effects, and to address ellipticity- 
gradient bias within the Deimos method. A C++ implementation of 
the method described here can be requested from the authors. 
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